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Quantization of constraint systems within the Weyl-Wigner-Groenewold-Moyal framework is dis- 
cussed. Constraint dynamics of classical and quantum systems is reformulated using the skew- 
gradient projection formalism. The quantum deformation of the Dirac bracket is generalized to 
match smoothly the classical Dirac bracket in and outside of the constraint submanifold in the limit 
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I. INTRODUCTION 

Gauge symmetries provide mathematical basis for 
known fundamental interactions. Within the generalized 
Hamiltonian framework 1] , gauge theories correspond to 
first-class constraints systems. Upon gauge fixing, these 
systems convert to second-class constraint systems. The 
operator quantization schemes for constraints systems 
have been developed by Dirac 0. The path integral 
quantization has also been developed and found to be es- 
pecially effective for gauge theories (for reviews see 0,0)- 

Besides conventional operator formulation of quantum 
mechanics and the path integral method, the popular 
approach to quantization of classical systems is based on 
the Groenewold star-product formalism^. It takes the 
origin from the Weyl's association rule [5( between oper- 
ators in the Hilbert space and functions in phase space 
and the Wigner function 0. The star-product formal- 
ism is known also under the names of the deformation 
quantization and the Moyal quantization 0] ■ 

The skew-symmetric part of the star-product, named 
the Moyal bracket, governs the evolution of quantum sys- 
tems in phase space, just like the Poisson bracket governs 
the evolution of classical unconstrained systems and the 
Dirac bracket governs the evolution of classical constraint 
systems. The Moyal bracket represents the quantum de- 
formation of the Poisson bracket. The quantum deforma- 
tion of the Dirac bracket has been constructed recently 

a 

The outline of the paper is as follows: In the next 
Sect., we give a pedagogical introduction to the Weyl's 
association rule using the elegant method developed by 
Stratonovich |Toj and give an introduction to the star- 
product formalism. More details on this subject can be 
found in articles [H El El El El 13 . 

The phase-space functions and the Dirac bracket do 
not make any physical sense outside of constraint sub- 
manifolds. In Ref. f| we constructed the quantum de- 
formation of the Dirac bracket on the constraint subman- 



ifold, sufficient for the purpose of generating time evolu- 
tion of quantum constraint systems. It would, however, 
be interesting from the abstract point of view to have 
a quantum-mechanical extension of the Dirac bracket 
which matches smoothly at ft — ► with the classical Dirac 
bracket outside of the constraint submanifold also. 

This problem is addressed and solved in Sects. Ill 
and IV. In Sect. Ill, we reformulate the classical con- 
straint dynamics using projection formalism and present 
the classical Dirac bracket of functions in terms of the 
Poisson bracket of functions projected onto constraint 
submanifold. Sect. IV gives the quantum-mechanical 
generalization of the method proposed. Sects. III-D and 
IV-B,C contain new results, the others is a pedagogical 
exposition of earlier works (mainly [9|). 

In Conclusion, wc summarize results. 



II. WEYL'S ASSOCIATION RULE AND THE 
STAR-PRODUCT 

Systems with n degrees of freedom are described 
by 2n canonical coordinates and momenta £ l = 
(q 1 , . . . , q n , pi , . . . , p n ) . These variables parameterize phase 
space T»R™ defined as the cotangent bundle of n- 
dimensional configuration space W 1 . Canonical variables 
satisfy the Poisson bracket relations 



= -i kl - 

The skew-symmetric matrix I kl has the form 

-E n 

E n 



(11.1) 



(II.2) 
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where E n is the n x n identity matrix and imparts to 
T*R ra a skew-symmetric bilinear form. The phase space 
acquires thereby structure of symplectic space. The dis- 
tance between two points in phase space is not defined. 
One can measure, however, areas stretched on any two 
vectors £ fe and as A = hit; k ( l where I kt — —I kl so 
that I k il lm = 

Principal similarities and distinctions between Eu- 
clidean and symplectic spaces are cataloguized in Table 
1. For skew-gradients of functions, short notation Idf(£) 
is used. 
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TABLE I: Comparison of properties of Euclidean and sym- 
plectic spaces 



[h i i /-• 1 1 /-J g-\ n y% T"\ O f~* /~\ 

JZj LlCllUfcJtlll oJJdCfc; 




x,y G R n 


f, c e R 2n 


Metric structure 


Symplectic structure 


<7ij = <?js 

„ n jk cfc 

ftjff — °i 


lij ~ —Iji 
I ij 1 — °i 


Scalar product 


Skew — scalar product 


(x,y) = Qij^y 3 


(00 =iii?C j 


Distance 


Area 


L = v( x -y, x -y) 


-4 = (00 




SlfPW — PTflHlPTTi" 


q(vfT = H <>1 <''■'■ 


(7d/) i = -I iJ df/d£ j 
= {CJ} 


Scalar product 


Poisson bracket 


of gradients of/ and g 


of / and g 


(vf, V<?) 


(Idf,Idg) = {f,g} 


Orthogonality 


Skew — orthogonality 


gijxY = 


ine< j = o 



In quantum mechanics, canonical variables are as- 
sociated to operators of canonical coordinates and mo- 
menta f = (q 1 , q", pi, p n ) acting in the Hilbert 
space, which obey the commutation relations 



[f,f] = -ihl kl . 



(11.3) 



The Weyl's association rule extends the correspondence 
<-> f to phase-space functions /(0 S C°°(T*ffi") and 
operators f 6 Op(L 2 (W 1 )). It can be illustrated as fol- 
lows: 



e TM n 

{0,0} 
/(0 e C°°(T*R n ) 



f € Op{L 2 (R n )) 
feOp(L 2 (K n )) 



The set of operators f acting in the Hilbert space is 
closed under multiplication of operators by c-numbers 
and summation of operators. Such a set constitutes vec- 
tor space: 



c x /(0 
/(0 + <?(0 

/(0*<?(0 



cf 

f+J 
ffl 



vector 
space 



algebra 



Elements of basis of such a vector space can be labelled 
by canonical variables 0- The commonly used Weyl's 
basis looks like 

35(0 = (27rft)"5 2 "(e-r) 



(IL4) 



The objects 33(0 satisfy relations 

33(0+ = 
Tr[B(0] = 
d 2n t 



(2tt?j) 

Tr[33(£)03(O] = (27rfi)"c5 2n (£ 

®(Oexp(-yP«0®(O 
= (27rfr)"<S 2 "(£- 0)03(0). 



1, 

1, 

f, 

a 



Here, 



<p = _rkl^_ ^_ 

is the so-called Poisson operator. 

The Weyl's association rule for a function /(0 and an 
operator f has the form 0] 



/(0 



I>[93(Of], 
d 2n £ 
(2ttS)' 



/(0®(0- 



In particular, 



Tr[?8(0? 
d 2n £ 



(2tt/i) 



(II.5) 
(II.6) 

(II.7) 
(II.8) 



The function /(0 can be treated as the coordinate of f 
in the basis 33(0, while the right side of Ea. (|II.5|l can be 
interpreted as the scalar product of 33(0 an d f- 

Alternative operator bases and their relations are dis- 
cussed in Refs. 0,^3- One can make, in particular, 
operator transforms on 33(0 an d c- number transforms 
on 0- Ambiguities in the choice of operator basis are 
connected to ambiguities in quantization of classical sys- 
tems, better known as "operator ordering problem". 

The set of operators is closed under multiplication of 
operators. The vector space of operators is endowed 
thereby with an associative algebra structure. Given two 
functions /(£) = Tr [03 (Of] and g(£) = 7>[33(Og], one 
can construct a third function 



/(0*s(0 = rrp8(0fe] 



(II.9) 



This operation is called star-product. It has been in- 
troduced by Groenewold £|. The explicit form of the 
star-product is as follows: 

/(0*ff(0 = /(0exp(jP)3(0. (H.10) 



where V = Vm- 
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The star-product splits into symmetric and skew- 
symmetric parts 



III. CLASSICAL CONSTRAINT SYSTEMS IN 
PHASE SPACE 



f*g = f°9 + ^-fAg. 



(11.11) 



The skew-symmetric part / A g is known under the name 
of Moyal bracket. It is essentially unique [l?J. It governs 
quantum evolution in phase space and endows the set of 
functions with the Poisson algebra structure: 

physical observables 
t 

functions in phase space 



Poisson algebra 



(11.12) 



f + 9, cx f, f*g,fAg 



vector space 



algebra 



The average values of a physical observable described 
by function /(£) are calculated in terms of the Wigner 
function 



W(0=Tr[«B(Ot]. 



It is normalized to unity 
d 2n £, 



w(t) = i. 



(11.13) 



(11.14) 



Iff- 
then 



/(£) and r <-» W(£) where r is the density matrix, 



Tr[fx 



d 2n £ 



(11.15) 



Under the sign of integral, the star-product can be re- 
placed with the pointwise product [l(j ■ 

Real functions in phase space stand for physical ob- 
servables, which constitute in turn the Poisson algebra. If 
the associative product f-kg does not commute, its skew- 
symmetric part gives automatically the skew-symmetric 
product which satisfies the Leibniz' law 



/ A(g*h) = (/ Ag)*h + g*(f Ah). 



This equation is valid separately for symmetric and skew- 
symmetric parts of the star-product. In the last case, 
Eq. _ 



(|II.16|) provides the Jacobi identity. The validity of 
the Leibniz' law allows to link the Moyal bracket with 
time derivative of functions and build up thereby an evo- 
lution equation for functions in phase space. 

In classical limit, the Moyal bracket turns to the Pois- 
son bracket: 

Hm/Aff = {f,g}, 



Second-class constraints Ga{Q = with a = 1, ...,2m 
and m < n have the Poisson bracket relations which form 
a non-degenerate 2m x 2m matrix 



(III.1) 



If this condition is not fulfilled, it would mean that gauge 
degrees of freedom appear in the system. After impos- 
ing gauge-fixing conditions, we could arrive at inequality 
Hill. 10 . Alternatively, breaking condition (fill. L|> would 
mean that constraint functions are dependent. After re- 
moving redundant constraints, we arrive at inequality 



Constraint functions are equivalent if they describe the 
same constraint submanifold. Within this class one can 
make transformations without changing dynamics. 



A. Symplectic basis for constraint functions 

For arbitrary point £ of constraint submanifold T* = 
{£, : G a (0 — 0}j there is a neighborhood where one may 
find equivalent constraint functions in terms of which the 
Poisson bracket relations look like 



where 



E n 

—E m 



(III.2) 



(III.3) 



Here, E m is the identity m x m matrix, X a fcZf, c = — 6 ac . 

The global existence of symplectic basis (|III.2() is an 
opened question in general case. The basis (|III.2|) always 
exists locally, i.e., in a finite neighborhood of any point 
of the constraint submanifold. This is sufficient for needs 
of perturbation theory. The formalism presented in this 
section can therefore to be used to formulate evolution 
problem of any second-class constraints system in phase 
space in the sense of the perturbation theory. 
The existence of the local symplectic basis 



(jrrpii is 

on the line with the Darboux's theorem (see, e.g., [lSjl 
which states that in symplectic space around any point 

(11.16) 

£ there exists coordinate system in such that £ £ 
where symplectic structure takes the standard canonical 
form. Symplectic spaces can be covered by such coordi- 
nate systems. 

This is in contrast to Riemannian geometry where 
metric tensor at any given point x can always be made 
Minkowskian, but in any neighborhood of x the variance 
of the Riemannian metric with the Minkowskian metric 
is, in general, ~ Ax 2 . Physically, by passing to iner- 
tial coordinate frame one can remove gravitation fields 
at any given point, but not in an entire neighborhood 
of that point. The Darboux's theorem states, reversely, 
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that the symplectic structure can be made to take the 
standard canonical form in an entire neighborhood 
of any point £. In Riemannian spaces, locally means at 
some given point. In symplectic spaces, locally means at 
some given point and in an entire neighborhood of that 
point. 

Locally, all symplectic spaces are indistinguishable. 
Conditionally, one can say that any surface in symplectic 
space, including any constraint surface, is a plane. 

In the view of this marked dissimilarity, the validity of 
Eas. ljII.lfl in a finite domain looks indispensable. 



B. Skew-gradient projection 

The concept of skew-gradient projection £ s (£) of 
canonical variables £ onto constraint submanifold plays 
very important role in the Moyal quantization of con- 
straint systems. Geometrically, skew-gradient projec- 
tion acts along phase flows IdQ a (J;) generated by con- 
straint functions. These flows are commutative in virtue 
of Eas. ljIII.2fl : Using Eas. ljIII.2f) and the Jacobi identity, 
one gets {G a ,{g b J}} = {G b ,{G a J}} for any function 
/, so the point of intersection with T* is unique. Skew- 
gradient projections are investigated in Refs. |l9j | and 
independently in Refs. 0, ■ 



in the power series of Q a , one gets 

^(0 = E^{."{tt,e ai },a aa },...e°*} 

fe=0 

><G ai Ga 2 -Ga k - (m.6) 

Similar projection can be made for function /(£): 

oo 1 

k=0 

xGarGaa-Ga*- (HI.7) 

It satisfies 

/ s (0 = f(UO)- (ni.8) 

Constraint functions are in involution with projected 
function: 

{/.(0.3»(0} = o. (ni.9) 

Consequently, / s (£) does not vary along IdQ a {^), since 

{f(o,9(0} = ^W-(id9(oy- 




FIG. 1: Schematic presentation of skew-gradient projection 
onto constraint submanifold along commuting phase flows 
generated by constraint functions. 

To construct skew-gradient projections, we start from 
equations 

{&(0.0«(O} = o (iii.4) 

which say that point £ s (£) 6 T* is left invariant by phase 
flows generated by G a (Q- Using symplectic basis (|III.2|) 
for the constraints and expanding 

UO = Z + xa Ga + \x ab g a g b + ... (ins) 



Applying Eas. (jlll.7fl and (|111.8|l to constraint functions 
Ga(Q, one concludes that the point £ s (£) belongs to the 
constraint submanifold 

Ga(UO)=0- (IILIO) 

The constraint submanifold can therefore be described 
equivalently as T* = {£ s (£) : £ E T*M n }. 

An average of function /(£) is calculated using the 
probability density distribution p(£) and the Liouville 
measure restricted to the constraint submanifold |2l| : 

^(2n) m l[8(Ga(li))mP(0- (IILll) 

On the constraint submanifold £ s (£) = £> so /(£) and 
p(£) can be replaced with /«(£) and p s (£). 

There exist therefore equivalence classes of functions 
in phase space: 

/(0~ff(0~/.(0 = s.(0- (ni.12) 

The symbol ~ means that functions are equal in the 
weak sense, /(£) rs <?(£), i.e., on the constraint submani- 
fold. We shall see that symbols ~ and « acquire distinct 
meaning upon quantization. Note that /(£) ~ /s(£)- 
Eas. HIII.8|l and (|III.10|) imply Q a ~ 0. Constraint func- 
tions belong to an equivalence class containing zero. 
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C. Dirac bracket in terms of Poisson bracket on 
constraint submanifold 

Given hamiltonian function Ti., the evolution of func- 
tion / is described using the Dirac bracket 



2m 



d_ 

dt 



f = {f,n} D . 



(III. 13) 



In the symplectic basis (|III.2|) . the Dirac bracket looks 
like 



Lc a ldG a {$)\ I IdH($ 

a=l 




T* = {£:G(Q = 0} 



{/, 9}d = {/, g} + {/, Q a }{Q a , g}- (ni.14) 

On the constraint submanifold, one has 

{f,g}n = {f,g s } = {f s ,g} = {f s ,g s }. (HI.15) 

Calculation of the Dirac bracket can be replaced there- 
fore with calculation of the Poisson bracket for functions 
projected onto the constraint submanifold. 

Two functions are equivalent provided they coincide on 
the constraint submanifold. The hamiltonian functions 
determine the evolution of systems and play thereby spe- 
cial role. Two hamiltonian functions are equivalent if 
they generate within T* phase flows whose projections 
onto the tangent plane of the constraint submanifold are 
identical. One may suppose that the equivalence relation 
for functions, defined above, does not apply to hamilto- 
nian functions, since skew-gradients of hamiltonian func- 
tions enter the problem cither. This is not the case, how- 
ever. The components of the hamiltonian phase flow, 
which belong to a subspace spanned at T* by phase flows 
of the constraint functions, do not affect dynamics and 
could be different, whereas the skew-gradient projection 
(|III.7|I does not modify components of skew- gradients of 
functions, tangent to constraint submanifold. We illus- 
trate it schematically on Fig. 2. The geometrical sense 
of the Dirac bracket reduces to dropping the component 
of the hamiltonian phase flow which does not belong to 
tangent plane of the constraint submanifold. Equiva- 
lently, those components can be made to vanish with 
the help of the skew-gradient projection. TL and TL S are 
thereby dynamically equivalent, so Eci. (|III.12fl character- 
izes an equivalence class for the hamiltonian functions ei- 
ther. Among functions of this class, TL S is the one whose 
phase flow is skew-orthogonal to phase flows of the con- 
straint functions, i.e., {G a ,T~Ls} = (IdQ a , IdTC s ) = 0. 

Replacing TL with TL S , one can rewrite the evolution 
equation in terms of the Poisson bracket (cf. Ea. (|III.13|l h 

^f = {f,n s }. (ni.16) 

The evolution does not mix up the equivalence classes. 

The physical observables in second-class constraints 
systems are associated with the equivalence classes of 
real functions in the unconstrained phase space. The 



FIG. 2: Schematic presentation of phase flows IdTL(£.) and 
IdTL s (£) generated by hamiltonian function H(£) and pro- 
jected hamiltonian function H. s (0 at point £ of constraint 
submanifold F*. The phase flow IdTL s (£,) belongs to the tan- 
gent plane of F*. The hamiltonian phase flow IdTL((,) admits 
decomposition IdH(0 = E„=i c a Idg a {S) +IdH s (£). Within 
the constraint submanifold (i.e. £ £ T* and £ + d£ 6 T*) 
one has dg a (£) = and therefore = dCdG a {£)/dC = 
(Idg a (£),dO- The first term E^i^WK) is therefore 
skew-orthogonal to any vector d£ of the tangent plane. 

equivalence classes constitute a vector space O equipped 
with two multiplication operations, the associative point- 
wise product and the skew- symmetric Dirac bracket {, }u, 
which confer O a Poisson algebra structure. 

Instead of working with equivalence classes of functions 
£f, one can work with their representatives f s defined 
uniquely by the skew-gradient projection. The one-to- 
one mapping £f <-> f s induces a Poisson algebra structure 
on the set of projected functions. The sum £ f + £ g con- 
verts to f s + g s , the associative product £f£ g converts 
to the pointwise product f s g s , while the Dirac bracket 
becomes the Poisson bracket: 

!/-•;/- I n !/-'/-!. (HI.17) 

These operations satisfy the Leibniz' law and the Jacobi 
identity and, since (f s + g s ) s = f s + g s , (f s g s )s = f s g s , 
and {f s ,h s } s — {f s ,h s }, keep the set of projected func- 
tions closed. 



D. Dirac bracket in terms of Poisson bracket on 
and outside of constraint submanifold 

Outside of the constraint submanifold functions do not 
make any physical sense. It is sufficient thus to work 
with the Dirac bracket on the constraint submanifold. 
The evolution problem in such a case can consistently be 
formulated in terms of the Poisson bracket for functions 
projected onto the constraint submanifold. 

The Dirac bracket is, however, well defined in the 
whole phase space. Redefinition of constraint functions 
by shifts G a {£.) — > Qa (£) + constant leaves the Dirac 
bracket unchanged, because it depends on derivatives of 
constraint functions only. It is not the case for the Pois- 
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son bracket applied to projected functions. This is why 
Eci. (|111.15|) is valid on constraint submanifold only. 

One can modify projection formalism to fit the above- 
mentioned property of the Dirac bracket. Suppose we 
wish to find the Dirac bracket of functions /(£) and 
at a point £ = £ outside of the constraint submanifold. 
The intersection of level sets {£ : Ga(() = Ga(0} can 
be considered as new constraint submanifold defined by 
constraint functions 

Projected functions depend thereby on both £ and £: 
MO = £^{...{{/(C),A(? Ql },A^},-A^} 

k=Q 

x Ag ai AG a2 ...AG ak (111.18) 

and similarly for <?(£)• The Poisson brackets are calcu- 
lated with respect to £ while £ is a parameter. The ap- 
propriate extension looks like 

{/(£), <K£)}d = {/(C),ss(C)}k=€ 
= {/s(C),<?(C)}lc=« 
= {/s(C),5s(C)}k=e- (in.19) 

In Eg. MI. 15(1 all four terms are pairwise distinct func- 
tions in the whole phase space. These functions coincide 
on the constraint submanifold only. In Ea. (|III.19jl all 
four terms coincide in the whole phase space. If £ G T* , 
we reproduce the result (|III.15|I derived earlier. 

IV. QUANTUM CONSTRAINT SYSTEMS IN 
PHASE SPACE 

Scheme presented in the previous Sect. 3 is suitable 
to approach description of quantum constraint systems 
in phase space. We give final results and refer to 9] for 
intermediate steps. 

We remind that classical hamiltonian function 7i(£) 
and constraint functions Ga(£,) are distinct in general 
from their quantum analogues H(£) and G a (£). These 
dissimilarities are connected to the usual ambiguities in 
quantization of classical systems, being not specific for 
the problem we are interested in. It is required only 

limff(0 = W(0, 

limG a (0 = G a (0- 

In what follows r* = {£ : G Q (£) = 0}. 

A. Quantum deformation of the Dirac bracket on 
constraint submanifold 

The quantum constraint functions G a (£) satisfy 

G a ($AG b (0=l ab - (IV. 1) 



In classical limit, G a (£) turn to G a (0- 

The quantum-mechanical version of the skew-gradient 
projections is defined with the use of the Moyal bracket 

6(0 A G a (0=0- (IV.2) 

The projected canonical variables have the form 

oo 

&® = J2^(- MAGai)AGa2) - AGak) 

k=0 

oG ai oG a2 ...oG ak . (IV.3) 
The quantum analogue of Eq. 1(111.7(1 is 

oo 1 

fc=0 

oG ai oG a2 ...oG afc . (IV.4) 

The function /*(£) obeys equation 

A(OAG tt (0 = 0. (IV.5) 

The evolution equation which is the analogue of 
Eg. ((III. 16(1 takes the form 

^/(0 = /(0AH t (0 (IV.6) 

where H t (£,) is the hamiltonian function projected onto 
the constraint submanifold as prescribed by Eg. ((IV. 4(1 . 
Taking projection of Eo. ((IV.6() we get evolution equation 
in the closed form for projected functions: 

^/ t (f) = /tK)AfftK) (iv.7) 

The quantum deformation of the Dirac bracket rep- 
resents the Moyal bracket for two functions projected 
quantum- mechanically onto the constraint submanifold. 

The formal structure of the dynamical quantum system 
is described by the scheme ((11.121 with the word "func- 
tions" replaced by the phrase "projected functions" and 
/ and g replaced by ft and <?t, respectively. The star- 
product is an associative operation, whereas the Moyal 
bracket for projected functions satisfies the Leibniz' law 
and, respectively, the Jacobi identity. 

Projected functions in phase space are objects associ- 
ated to quantum observables. Functions which have the 
same projections are physically equivalent. We can unify 
such functions into equivalence classes. The star-product 
and the Moyal bracket for projected functions generate 
for equivalence classes a Poisson algebra structure ac- 
cordingly. 

The bracket f t A g t constructed in |9| gives the deforma- 
tion of the Dirac bracket on T* . What about the whole 
phase space? 



7 



TABLE II: Brackets which govern evolution in phase space 
of functions (second column) and projected functions (third 
column) of classical systems (first row) and quantum systems 
(second row) . The right upper corner shows the Dirac bracket 
expressed in terms of the Poisson bracket of functions pro- 
jected onto the constraint submanifold. The left upper cor- 
ner is the Poisson bracket. The left lower corner is the Moyal 
bracket, which represents the quantum deformation of the 
Poisson bracket. The operation ft A gt is the quantum defor- 
mation of the Dirac bracket. 



Systems: 


unconstrained 


constrained 


classical 


lf,g\ 


{fs,9s} 


quantum 




ft Afft 



B. Quantum deformation of the Dirac bracket on 
and outside of constraint submanifold 

One can generalize the operation f t A g t to match in 
classical limit the Dirac bracket outside of the constraint 
submanifold. We can proceed like in the classical case by 
writing projected functions in the form 

MO = ]T M (...((/(C)AAG ai )AAG a2 )A 

fc=0 

... AG ak ) o AG ai o AG Q2 ... o AG afc (IV.8) 

where 

AG o (C)=G o (0-G o (£). 

The Moyal brackets and the o-products entering this 
equation are calculated with respect to (. The desired 
extension looks like 

/r(C)A0r(Olc=e- ( IV -9) 

It is assumed that the constraint functions G a (£) satisfy 
the bracket relations HIV. ID at £ ^ T*. Expression HIV.9|) 
is valid on and outside of the constraint submanifold. If 
£ 6 r*, we reproduce operation /*(£) A <?t(£) announced 
earlier. 



C. Completeness of the set of projected operators 
of canonical coordinates and momenta 

The set of operators f is known to be complete, so 
that any operator f can be represented as a symmetrized 
(probably infinite) weighted sum of products of opera- 
tors f. In the sense of the Taylor expansion, one can 
write f = /(y). The one-to-one correspondence between 
operators f € Op(L 2 (R n )) and functions in phase space 
/(£), based on the Taylor expansion, is equivalent to the 
Weyl's association rule. 



The similar completeness condition holds for projected 
operators of canonical variables which are inverse 
Weyl's transforms of Apparently, any operator 

f acting in the Hilbert space can be represented as an 
operator function (p(& a ,f t ). Applying projection to the 
symmetrized product of k constraint operators Q5 a , which 
are inverse Weyl's transforms of G a (£), one gets a series 
like l — k+ ±k{k - 1) + . . . = (1 - l) fc = 0, and so 

(<g(<*i<ga 2 ...0»*)) t = o. (IV. 10) 

The Taylor series of <p(<S a , f t ) generates thereby vanishing 
terms involving (3 a . We thus obtain 

M0 a , y j)) t = ^(o, t j). (rv.ii) 

Respectively, any function projected quantum- 
mechanically onto the constraint submanifold can 
be represented in the form 

MO = y(*6(0)- (iv. 12) 



One can pass to classical limit to get Ea. (IIII.8ll . Con- 
structing <p(£) from /(£) is a non-trivial task equiva- 
lent to solving constraints. The operator counterpart of 
Ea. (ITVT2l) . 

U = <p(ft), (IV.13) 

demonstrates the completeness of projected set of oper- 
ators of canonical coordinates and momenta. Accord- 
ingly, Eg. HIV. 12(1 shows completeness of the set of £|(£) 
in description of projected functions. It is worthwhile to 
notice that Eq. i|lV.10(l does not extend to antisymmetric 
products of <3 a as one sees from [6 a ,6 b ] t = {-l ab ) t = 
-l ah ^ [0?,©?] = where condition 0° = is taken 
into account. 



V. CONCLUSION 

We made short introduction to the Weyl's association 
rule and the Groenewold star-product technique for un- 
constrained and constraint systems. The attention was 
focused to the evolution problem. 

A generalization of the quantum deformation of the 
Dirac bracket is constructed to match smoothly classical 
Dirac bracket in the whole phase space at K — > 0. 

The use of skew-gradient projection formalism allows 
to treat unconstrained and constraint systems essentially 
on the same footing. Projections of solutions of quan- 
tum evolution equations onto the constraint submanifold 
comprise the entire information on quantum dynamics of 
constraint systems. 
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